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We extend the dimensional deconstruction by utilizing the knowledge of graph theory. In the 
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solutions in our model. We show some examples of the graphs admitting the vortex solutions of 
simple structure as the Bogomolnyi solution. 

PACS numbers: 02. 10. Ox, ll.10.Lm, 11.27.+d, 11.30.Qc 



*Elcctronic address: kan@yamaguchi-jc.ac.jp 
^Electronic address: m004wa@yamaguchi-u.ac.jp 
^Electronic address: shiraish@sci.yamaguchi-u.ac.jp 



1 



I. INTRODUCTION 



0,3- 



Most of these 



Recently, 'Higgsless theories' gerly studied by many authors 

models are derived from or related with the method of the dimensional deconstruction 
(DD) [3], which leads to the breakdown of electroweak symmetry. 

The typical structure of DD is shown diagrammatically in FIG. [TJ This model incorpo- 
rates the [SU(2)] N+1 ® U(l) gauge group and N + 1 nonlinear-sigma-model fields. If N is 
equal to one, the number of the site is three in FIG. HJ The three-site Higgsless model 
is in this category. In the generic scenario, the [SU(2)] N+l ® U(l) gauge group is broken to 



FIG. 1: A moose diagram. There are N + 1 SU(2) gauge fields and a U(l) gauge field. Each gauge 
field exists on each site represented by a small circle. The coupling constant of the gauge fields 
A 1 is gi (i = 0, 1, 2, • • ■ , N), and the coupling constant of the U(l) gauge field is gN+i- The 
vacuum expectation value of the scalar fields Ej is fi (i = 1, 2, • • ■ , N, N + 1). 

The moose diagram like FIG. [1] naturally leads to the Lagrangian of the model. This 
moose diagram indicates a relation between gauge fields and scalar fields. We will generalize 
this relation in the context of graph theory. We can express the relation between gauge 
fields and scalar fields in a graph, which is just a complex moose. We wish to call this 
theory based on a graph as 'graph dimensional deconstruction' (GDD). The idea of GDD 
has already been published as Ref. 4|. 

In the present work, we propose another idea of using superfields to introduce supersym- 
metry (SUSY) into the model. We assign vector superfields to vertices and chiral superfields 
to edges of a graph. This is another extension of the DD. 

In the beginning, both DD and SUSY are to provide the mechanism of solving the gauge 
hierarchy problem. The motivations of including SUSY are, nevertheless, claimed as follows. 
First of all, we should think that every field theory has SUSY at very high energy, because 
the correct or controlled UV behaviors are believed, or because of superstring theory or 



M-theory. The second motivation comes from the necessity of more symmetries. Because 
DD and GDD are basically the mechanism of controlling the mass spectrum of field theory, 
we need more symmetry to determine the (self-)interaction of fields. Thus we consider the 



In this paper, we consider only the Abelian theory. For notation, please consult Ref. [4]. 

II. A REVIEW OF FIELD THEORY ON A GRAPH (OR GDD) 

A graph G(V, E) consists of a set of vertices V and a set of edges E. A vertex is connected 
with another one by an edge. We let the number of the vertices be p, p = j^V and the number 
of the edges be q, q = #E. In FIG. El we show the simplest graph with p = 2 and q = 1, 
constructed by two vertices and an edge. 



FIG. 2: The simplest graph, constructed by two vertices and an edge. A vertex Vi is identified 
by i, where i is a label for each vertex. In the same way, an edge is identified by i, where i is 
a label for each edges. The arrow means a direction of the edge. This edge is called an oriented 
edge. In terms of the oriented edge, the original vertex v\ is v\ = o{e\) and the terminal vertex vi 
is t>2 = t(ei). This oriented graph corresponds to the generalized moose diagram. 

We consider a simple Abelian theory. Abelian gauge fields reside at vertices and scalar 
fields reside at edges. The U(l) transformation is defined at each vertex. The Lagrangian 
density is 



super symmetric extension of the GDD model here. 



vi - o(e 1 ) v 2 - tie!) 





(i) 



where the covariant derivative is 



Ty>U e = (P + igAZ e) -igA% e) )U t 



(2) 



with \U e 



f 2 - 
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If we rewrite U e as U e = /e me , the real scalar fields a e act as the Stueckelberg fields. 
The number of physical massless scalar fields is q — p + 1, or the number of closed circuits 
involved in the graph, because p — 1 scalar degrees of freedom are absorbed by the to-be 
massive vector fields. If and only if the graph is tree (or absent from closed circuits), the 
scalar fields disappear from the physical spectrum. 

The (mass) 2 matrix of vector fields M\ is given by 2g 2 f 2 A, where the (p, p) matrix 

A = EE T , (3) 

n 

is called as the graph Laplacian and the (p, q) matrix E is the incidence matrix |16j] defined 

as 

1 if v = o(e) 

(E) ve = I -1 if v = t(e) ■ (4) 
otherwise 

Here v = o(e) means that the vertex v is the origin of the edge e and v = t(e) means that 
the vertex v is the terminus of the edge e. The (q, p) matrix E T is the transposed matrix of 
E. 

For more general cases, one might consider individual coupling constants for vertices as 

V»U e = (0" + ig t ( e) A^ e) - ig o{e) A^ e) )U e , (5) 
and \U e \ 2 = f 2 for each edge. In this case the mass matrix becomes 

M 2 A = 2GEF 2 E T G = 2(GEF)(GEE) T , (6) 
where the diagonal matrices G and F are given by 

\ Qv if v = v' f e if e = e' 

(G)w' ={ y , (F) ee , = \ Je , (7) 

otherwise otherwise 

respectively. 

To summarize this section: In the GDD model, the mass spectrum is given by eigenvalues 
of the graph Laplacian or the related matrix constructed from the incidence matrix of the 
graph. 
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III. THE USE OF THE STUECKELBERG SUPERFIELD 

Next we incorporate SUSY into the GDD model. We use superfields js] to this end. 
In this paper, we consider that vector superfields {V v } exist on vertices. We still impose 
the U(l) transformation on V v at each vertex as 

V V ^V V + i(K - A v ) , (8) 

where is a chiral superfield. Then the invariant superfield is defined as usual: 

W: = ~DDD a V v . (9) 

The kinetic term of the vector field can be created from this for each vertex. 

Further we introduce a chiral superfield S e at each edge. The superfield S e is assumed to 
be transformed as: 

S e -> S e - zA f ( e ) + iA (e) • (10) 

Then we can write the Stueckelberg term Q 

(V t(e) -V o{e) + S e + S e ) 2 , (11) 
and a gauge invariant term for the interaction with scalars 

C = E h { W v W a\o0 + ^WtL) + E 2 /e W« - V °to + Se + ^ e)\ m ■ (12) 

The bosonic part of the theory is found to be 

A = " E h F l^ ~ E ^Ke) ~ A o(e) + d^ef ~ \ E 2 /eV Pe f 
v£V yv egB eeE 

+ E + 2 E 2 /e l^el 2 + E 2 /e(A(e) " A*.))* , (13) 

v£V ® v egB eeB 

where the notation of component field is rather standard one and is gathered in Appendix [Al 
Eliminating the auxiliary fields F Se and rescaling p e , gauge fields and D v to have canonical 
kinetic terms we get 

f2 



ueV egB egB 

- E E^( eT )-^(£)-^^^ . (M) 

e,e'gBugV v&V I egB J 
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Now one can easily find the mass matrices for vectors and scalars: 

M\ = 2GEF 2 E T G = 2{GEF){GEF) T , M 2 = 2FE T G 2 EF = 2(GEF) T (GEF) , (15) 

where E is defined as (j3J) while G and F are given by (jTJ). Massless scalar fields are absent 
if and only if the graph is a tree graph. The mass spectrum of the scalar fields is the same 



as the one for the vector fields except for zero modes. 
The fermionic part of the theory is found to be 



veV ® v ee-B 
+ E 2 ^ U^Me) ~ Me)) + h.c] , (16) 



e&E 

and can be rescaled as 

[feXe(E T ) ev9v X v + h.C.] . (17) 

e&E v€V 

Here A„ and Xe are Weyl spinor fields contained in V v and S e , respectively. 
One will find the mass matrices for fermions after rescaling the fields: 

M\ = 2GEF 2 E T G = 2(GEF)(GEF) T , M 2 = 2FE T G 2 EF = 2(GEF) T (GEF). (18) 

Note that the fermions A and x form Dirac fields for massive modes. Also note that all field 
contents are neutral as well as free from interactions. 



IV. MULTI- VECTOR, MULTI-HIGGS MODEL 
A. general construction 

We will construct the model that the symmetry [77(l)] p is spontaneously broken to U(l). 
Therefore we will not use the Stueckelberg fields but the Higgs fields. 

As the model in the previous section, we consider vector superfields on vertices and 
suppose that U(l) transformation is defined at each vertex. Moreover in the present case, 
we introduce a 'bi-charged' scalar field S on each edge, which is transformed under two U(l) 



symmetries as 



it 



e -2*A t(e) Sg e 2*A o(e) _ ( 1Q ) 



Now we get the [?7(l)] p invariant supersymmetric multi-vector, multi-'Higgs' model on a 
graph governed by the following Lagrangian: 



- 2gJ2 Ce(V tie) -V o{e) )\ em , (20) 

e&V 

where we rescale the gauge coupling constant to be seen explicitly. The Fayet-Illiopoulos 
terms are chosen so that they are similar to those in the model of the previous section, when 
Ce ^ /e ■ 19] This paper will not go into the issue about anomaly and deal with only classical 
aspects of the model. 

The bosonic part of the Lagrangian reads 



a = -\J2 F ^ F " V+ IJ2 D "-J2( V ^( V ^ 

+ J2 F L F *e + <^(A(e) - D o{e) )a\a e - g £(A(e) - D o{e) )C e , (21) 



eG-E e£E e£V 

where the covariant derivative is 



V»a e = (d» + igA>t ie) -igA» ie) )a e . (22) 
By use of the incidence matrix of the graph, we rewrite the above Lagrangian as 

veV ueV e£E 

+ E^e^e-^E(^e-Ce)(^)e. (23) 



eeE eeE 

Substituting the equation of motion for the auxiliary fields 



F Se = and D v = gJ2(vl*e-(e)(E T ) v , (24) 



e€E 

into the bosonic Lagrangian, we obtain 



4 

vGV e£E 

9 - E (°l°e - Q{E T E) ee ,{ol,a e , - Q) . (25) 

e,e'€E 



Note that E T E is a (q, q) matrix. 
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Vl 




pA 



pB 



FIG. 3: P3: the path graph with three vertices. There are two substantially different graphs. They 
have the different incidence matrices. 

B. example: P3 



The structure of the model depends on the incidence matrix of the graph. For a simple 
example, let us consider the path graph with three vertices, P3. 

The incidence matrix depends on the orientation of edges. For instance, two cases can 
be considered as follows: [20J 



( 1 \ 



Alve 



V 



-1 1 
-1 



1 1 \ 



B ve 



\ 



-1 -1 

1 



where Ea is the incidence matrix of P% and Eb is the one of P 3 S . 

Interestingly, the following matrix is independent of the edge orientation: 



1 1 -1 \ 

-1 2 -1 
-1 1 



A 



(26) 



(27) 



This is known as the graph Laplacian. 
On the other hand, we find 



E t a Ea 




E b E b 



2 1 
1 2 



(28) 



Therefore the shape of the 'Higgs' potential in Eq. (|2"5|) depends on the edge orientation. 
Figure H] illustrates the contour plots of the potentials in Eq. (J251) for the graphs P 3 A and 



pB 
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0.2 0.4 0.6 0.8 1 1.2 1.4 



FIG. 4: Contour plots of scalar potentials for the models based on P£ (left) and on (right), 
respectively. In both plots, potentials are normalized by g 2 / 4 , the contour spacing is 0.1, and the 
horizontal axis indicates |cri|// while the vertical axis indicates \o<2\j f . 

C. mass matrices for bosonic and fermionic fields 

Individually different gauge coupling constants will also be considered. The consequence 
of such consideration forces the bosonic part of the Lagrangian to be 

v£V eG-E 

- \ E - Ce)(£ T W(£WKW - CeO , (29) 

e,e'£E uel 7 

with 

V*<r e = {&* + igt(e)A? {e) - ig o{e) A^ e) )a e . (30) 

Here we assume that all ( e are positive and a/CI — fe- Thus the vacuum expectation value 
for |er e | is f e and physical scalar fields should be considered as the linear combinations of 
\a e \ — f e . Each phase part of a to-be massive scalar field is eaten by a vector field through 
the Higgs mechanism. Then the (mass) 2 matrices My for vector fields and M| for scalar 
fields in this case are 

M v = 2GEF 2 E T G = 2(GEF) (GEF) T , M 2 S = 2FE T G 2 EF = 2(GEF) T (GEF) , (31) 

where the matrices that appeared in the above equations are the same as (jlj) and (j7|). 

Although the shape of the potential with respect to \a e \ depends on the orientation of 
edges in the graph, the mass spectrum of the scalar fields is the same as the one for the 
vector fields except for zero modes, similarly to the model in the previous section. 



The number of the moduli of the potential is q — p + 1 for a general graph. This is equal 



to the number of independent closed circuits in the graph, 
expectation values of a e are determined rigidly if all ( e are positive. 
The fermionic part of the Lagrangian is 



2l| For tree graphs, the vacuum 



veV eG-E 

+ tyj2j2(°eME T ) ev g v \ v - at^ e (E T ) ev g v X) , (32) 

where A„ and ip e are Weyl spinor fields contained in V v and £ e , respectively. The covariant 
derivative on ip e is defined as T>^ip e = (<9 M + igt{e)-^t(e) ~ ^o(e)^o( e ))V , e- Substituting the 
vacuum expectation values (a e ) = f e , we find 

Ml = 2(GEF)(GEF) T , M| = 2(GEF) T (GEF) . (33) 

Since SUSY is unbroken, the bosonic and fermionic spectra are the same. 

In this paper, we have considered models with unbroken SUSY. The model with 'partially' 
broken SUSY is interesting, for some ( e < 0. The present analysis will not go into such 
models. 



V. VORTEX SOLUTION 

It is well known that the vortex solution can be found in the Abelian-Higgs model [8fl. 

CI 

In many papers, the solution is used as a simple model for a cosmic string |9[. We consider 
the vortex-type solutions in our model described in the previous section. 

Although an academic interest in our toy model is an important motivation for the 
following study, we also think that topological configurations are a key ingredient in recent 
studies in theoretical physics. A possibility is expected that a similar model provides an 
example of a complicated brane/string system. In the present paper, anyway, we study only 
simple vortex in our theory and their generalizations and possible applications to particle 
physics and cosmology are left for future work. 

Moreover we will consider only tree graphs as the bases of models. 
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A. Bogomolnyi equation 



In the Abelian-Higgs model, the vortex solution is well known Moreover, it is 

Be- 



known [lo| that supersymmetric U(l) theory satisfies the Bogomolnyi condition 
cause our model is also supersymmetric, the Bogomolnyi condition can be found. The 
equations of motion can be reduced to the following two sets of equations: 



^' = ^%,^(^ e (|cr e | 2 -C e ) 



and 



where i, j denote two spatial directions and e 1 ^ is the antisymmetric tensor. 
These equations are the Bogomolnyi equations. 
The energy per unit length of a vortex string can be written as 



(34) 



(35) 



8 



d x 



d 2 x 



^E F ^' + E(W(2>v e ) 

+ ^ E Ed ^ 2 - (e)(E T ) ev g 2 v (E) vel (\a e ,\ 2 - ( e ,) 

e,e'£Ev€V 

- A E {Fl j ± g v (E) ve (K\ 2 ~ CO} 2 + \ E ± ^Vja e 
± I E E l^g v (E) ve (e ~ i E ^dMlV^ 



(36) 



(37) 



For a solution of finite energy density, T>i<j e is equal to zero at spatial infinity. If the 
asymptotic behavior of o e is expressed by the azimuthal angle <p and an integer n e , i.e. 



(7, 



/ £ee ine *', the condition tells {E T g v A 



i Je 



n, 



,di(p, and then j d 2 x {E T Sijg v F^ 



Aim e . Therefore the energy density becomes 

S = [d 2 x - E {F l v j ± e^g v (E) ve (We\ 2 - Ce)} 2 + \ E \ v ^e ± ie ij V j( T e 

J L v€V e&E 

±27rE l^elCe- 



(38) 



e&E 



We deal with the lowest bound for the energy density read from this result. The vortex so- 
lution satisfying the Bogomolnyi equation (134113 51) has the energy density 2n J2 e eE l n e|Ce-|22 
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B. Bogomolnyi vortices and SUSY 



It is well known that the SUSY is partially broken in the topological background fields. 
Here we briefly describe the pattern of SUSY breaking in our model. Notation is indebted 
in According to SUSY, the variations of the gauginos A e are 

5 e \ v = ieD v + a^F v ^e. (39) 

Using the Bogomolnyi equations (1341) . and assuming the vortex string lies in the third di- 
rection for simplicity, the above variations are rewritten as 

5 € \ v = TiF^(l±a 3 )e. (40) 

This means that the half of the SUSY at the vertex is broken in the presence of the central 
magnetic flux of the vortex. 

The variations of partners of a e are 

5 e ij e = iV2ea^V fl a e , (41) 

where I" 1 ^ = d^ipe + K(9^)t( e ) ~ (s^ote))^ - ^ ^ e vor t ex string lies in the third direction, 
this reduces when the Bogomolnyi equations fl35l) hold, 

6 £ i/j e = iV2e {a l V x il) e + a 2 V 2 ip e ] = iV2e (a 1 ± io 2 )V x ^ e . (42) 

We find again that the half of the SUSY at the edge is broken in the presence of the magnetic 
flux. 



C. construction of vortices: ansatz 

Next we examine how we can obtain the explicit solutions in our model. For simplicity, 
we consider a common gauge coupling constant g and a single constant / = y/C,. In other 
words, we consider the case that G = gl and F = fl (where / is the identity matrix). 
Although we cannot tell about most general solutions, we take ansatz for simple, physically 
admissible type of vortex solutions. [23|] We impose the axially symmetric ansatz 

a e = p e (r)e in ^, (43) 

K = Pv(r) , (44) 
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on Bogomolnyi equations. Here we express the radial coordinate as r and the azimuthal 
angle as (p. The integers n e are winding numbers. The detailed calculation is shown in the 
Appendix [Ql We get the following Bogomolnyi equations, 

4 = _(,(*■*)-») 

Pe r 

^ = -»£(£),, (p 2 -/ 2 l. («) 

where the prime (') denotes the derivative with respect to r. These equations are the special 
case of the Bogomolnyi equations. 

D. examples of vortex solutions 

We show some concrete examples for the vortex solution in our model. To have the vortex 
solution we restrict the graph structure, or equivalently, the incident matrix E. Here we 
also consider configurations with the least winding numbers for simplicity and for feasibility 
in physical systems. 

We consider here the cases with the single-centered exact solution similar to the nor- 
mal vortex. The asymptotic behavior of general cases can be obtained and is shown in 
Appendix [El 

1. Example 1: P2 

The simplest case has two vertices and an edge. This graph is P2 graph. We show the 
graph in FIG. [51 




FIG. 5: i-*2 : the path graph with two vertices. 
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In this case, the incidence matrix and its transposed matrix are 

(E) ve = r), (E T ) ev =(l-l)- (47) 




Then considering the Bogomolnyi equations 



^ = -,5>)„ e (p 2 -/ 2 ) e , (48) 



Pe 



(gE T P-n) ( 



(49) 



the first one becomes 



% = -9 (P 2 ~ f) , (50) 
^ = +9 (P 2 - f) ■ (51) 

Therefore it is necessary to find a set of unique equations that we suppose the relation 
Pi{f) = —P2{r). On the other hand, in the second equation we notice 



j2(E T ) ev p v = (i -ij y X J Pi = 2Pi • (52) 

So, we get the following equations 

^ = -9 (P 2 ~ f) , (53) 



r 

pf_ = 2gP 1 - n 
p r 



(54) 



These equations can be reduced to 

P' 



- (p 2 - 1) , (55) 



p P-n 



(56) 



p x 

if we rescale the variables so that P(x) = 2gP\(r), p(x) = p(r)//, x = \/2gfr, n = 1 and 
the prime (') is the derivative with respect to x. These equations are precisely same as the 
normal Bogomolnyi equations. The normal Bogomolnyi equations is referred in Appendix 
The energy per unit length of the straight string is given by 2nf 2 in this case. General- 
ization to the case with the winding number n > 1 is trivial. 

14 



2. Example 2: P 3 



We consider the P 3 graph, the three-vertex path graph. In this graph, we consider two 
patterns of the direction of the edges. We show these in FIG. [HI 




FIG. 6: The graph P^ has edges of the same direction while P% has the edges of the different 
direction. 

The condition to reduce the Bogomolnyi equations in these cases to the normal ones 
( 15511561) with pi = P2 and n\ = ri2 = 1 are P\{r) = — Ps(r) and Pzir) = in the case with 
P% while P\{r) = P 3 (r) and P2(r) = —2P\(r) in the case with P 3 B . The necessary scaling is 
that P(x) = gP\{r) and x = gfr in the case with P 3 A while P{x) = 3gP\(r) and x = \^3gfr 
in the case with P 3 . The energy density takes the same value 27r/ 2 (l + 1) = 47r/ 2 in both 
cases. 

3. Exam-pie 3: 




FIG. 7: The star graphs, K^ N and K 
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We consider another tree graph, the star graph Ki^. In the star graph, vn+i is adjacent 
to all the other vertices and no extra edge exists. We recognize two types of edges. One is 
the edge whose origin is Ujv+i, another edge is one whose terminus is Ujy+i- We call the edge 
of the first type is e D , the one of the second type is e t . 

We heuristically find the cases that we get the vortex solution similar to the normal one 
with pi = pi = ■ ■ ■ = pn = pn+i- Here two cases are shown where the number of edges 
belonging to two types are 

Kt, N : #e = #e t = N/2, (57) 

Kin '■ # e o — N an d #e t = , or vice versa , (58) 

where, off course, N is considered to be even in the case A. The graphs of two types are 
shown in FIG. [7J 

The incidence matrix of K^ N (where N is even) is (N + 1, N) matrix given by 

( -\ ■■• \ 
1 ■■■ 
0-1 ■■■ 



(E. 



Ajve 



while the incidence matrix of K^ N is 



■■ 

V i -l i •• 

( -I •• 

-1 •■ 

-1 ■• 



(59) 



BJve 




V 1 





1 





1 



• - 1 / 

• o \ 

■ 

■ 

■ -1 

• 1 J 



(60) 



We found these patterns by extending the analysis of getting the vortex solution in the 
case with P 3 graph shown previously, because Kx 2 is the same as P3. 

In the first case flBTl) . we have vortex solutions if P2£-i( r ) = ~P2m( r ) (£,m are positive 
integers and £, m < y) and Pn+i = 0. In the second case (1581) . we have the solutions if 
Pi(r) = Piir) = ■ ■ ■ = P/v(r) and P/v + i(r) = — A r P i (r). In both cases the energy density is 
found to be 2TrNf 2 if all the winding numbers are unity. 
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4- Inclusion of 'no winding scalar edge' 



In the previous two examples, all 'Higgs' scalars have nonzero winding number. Con- 
versely we consider that there is an edge where the assigned scalar has no winding number, 
thus p e = / at the edge. We use the dashed line to express such an edge, as in FIG. |HJ 



FIG. 8: This dashed line means that p e = f on this edge, 'no winding scalar edge'. 



For a constant p e , P D ( e )(r) = P t ( e ){r) holds everywhere. \24i\ Suppose that one have already 
constructed the vortex solution in a certain model with specific graph structure. The one 
might duplicate the solution and the graph. One may connect the identical vertices of the 
original and copy of the graph by 'no winding scalar edge'. The number of such connection is 
arbitrary. This method can be applied to the case with two different models and solutions, if 
one finds the same functional form of P v (r) in each model. Of course more than two vertices 
can be connected if P v is common at all vertices. 

5. Example 4 : Pa 





FIG. 9: i-4 graph consists of two P2 and an edge. 

We consider the P4 graph. The graph P 4 has two Pi as subgraphs and is shown in FIG. 
We do not show the direction of the edge in this graph. This graph has a left-right symmetry 
with respect to the dashed edge. This symmetry is connected with the winding number of 
each vector fields. The vector fields at the both ends of the dashed line must be described 
by an identical function. For this reason, we should impose the left-right symmetry to the 
direction of edges. In the P4 case, we find two types of the edge orientation graph for 
admitting the normal vortex solutions, shown in FIG. [10] and FIG. [TTJ In the similar way, 
we consider the model based on P 2 £ with normal vortex solutions. 
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FIG. 10: P4 graph whose edge direction is left-right symmetric with respect to the dashed edge. 
Each of edge directions is outgoing with respect to the dashed edge. 




FIG. 11: P4 graph. Each of edge direction is incoming with respect to the dashed edge. 
6. Example 5: 

The graph Pq has three P2 as subgraphs. We study the model based on P$ and their 
standard solution in the above-mentioned way. In addition, P 6 has two P 3 as subgraphs. 




FIG. 12: Pg graph, which includes two P3 as subgraphs. 

Similarly to the case with P 4 , we can consider the P$ graph as two subgraphs connected by 
an edge. We exhibit the Pq graph in FIG. [12j We have the left-right symmetry with respect 
to the dashed edge also in this case. We classify four types of the graph in terms of the 
direction of the edges as in FIG. [131 In the similar way, we can consider the P^ graph, and 
associated models and solutions. 

7. Example 6 

We can connect two Ki >N graphs by the dashed edge as in FIG. [HI As this example, we 
can find the graph structure admitting the normal vortex solutions. 
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FIG. 13: There are four types of the Pq graph consisting of two P3. 



FIG. 14: The graph consisting of two 2£i,jv connected by the dashed edge. 

VI. CONCLUSION AND OUTLOOK 

We have generalized DD into GDD and introduced SUSY to GDD in the Abelian theory. 
A multi-Abelian-Higgs model has been studied as a further generalization. After getting 
the Bogomolnyi equations, we explicitly constructed vortex solutions of the normal type. 
To get the vortex solution, we restricted the graph structure to the special cases shown in 
the previous section. We showed some examples for the graph which has the normal vortex 
solution. 

We have left the following aspects of the multi-Abelian-Higgs models for future work. 
First, we discussed single-centered vortex in the present paper. The possibility of multi- 
vortex solution [l2|] is an important subject to study. Next, in this paper, we mainly con- 
sidered tree graphs. If we take general graph structures as the bases of multi-Abelian-Higgs 
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models, we have scalar potentials with (many) flat direction of the lowest energy. The ap- 
pearance of moduli is the feature of supersymmetric theories and the vortex solution in such 
a model is crucial for phenomenological models [ljj. At the same time, the quantum cor- 



rections might become essential. The generalization of the method in |14j will be useful to 
investigate the quantum effects about vortices. Finally, because our model contains several 
fields, the possibility of different types of topological defects, such as 'rings' [3], must be 
examined. 

We considered the Abelian gauge theory in GDD as well as multi-Higgs models. We are 
also interested in the non- Abelian theory because the three-site Higgsless model is based on 
the [S77(2)] 2 ® U(l) gauge theory. While we considered vortices in the Abelian gauge theory 
in this paper, on the other hand there exist monopoles in the non- Abelian gauge theory. As 
the future works, we wish to incorporate monopoles, superfields and GDD into non-Abelian 
theory as some toy models for the Higgsless model. 
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APPENDIX A: CONTENTS OF SUPERFIELDS 

In this Appendix, we collect the superfields and their component fields. See the refer- 
ence jf|. 

1. vector superfield 

v v = -QajdA* + ieeex v - mex v + -eeeeD v . (Ai) 



This satisfies 

1 
2 



V v 2 = --6666 A^A* , K 3 = • (A2) 
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2. chiral superfield (Stueckelberg superfield) 



1 -1 
S e = -(p e + ia e ) + 9xe + ida^9-{d ll p e + id tJL a e ) 

+B6F Se + l -eeea^ Xe + \eeee{ Pe + ia e ) , (A3) 



S e + S e = p e + 9 Xe + 0Xe - 9a»9d,a e + 99F Se + 00F] C 

+ l -999^d, X e + l -999^d,xe + \oOOO Pe . (A4) 



3. chiral superfield (Higgs superfield) 

S e = cr e + v^^e + i9a^9d^a e 

+99F Se + -^eeOa^d^e + -9999(a e ) . (A5) 
v 2 4 

APPENDIX B: THE EIGENVALUES OF MATRICES AB AND £A 

Let A be a (p, q) matrix and B be a (g, p) matrix. Then (p + q, p + q) matrices U and V 
are defined as 

u=( IrA )' < i ' ') < B1 > 

\B xl q j \ % p I q j 

where I p is the (p, p) identity matrix while qp is the (q,p) matrix all of which elements are 
zero. 

The products of two matrices are 

UV=\ Xlp %q ], VU=( Xlp ~ AB ° Pq I . (B2) 
xB xI q -BA J y B Xlg 

Because det UV = det VU, the eigenvalues of AB and BA are equal, except for zero eigen- 
values. 



APPENDIX C: THE NORMAL VORTEX IN ABELIAN-HIGGS MODEL 

The Ginzburg-Landau theory is used as a macroscopic theory of the superconductivity. 
That is nonrelativistic theory, and we know an Abelian-Higgs model as the relativistic version 
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of the Ginzburg-Landau theory. This model includes the normal vortex solution. In this 
paper we distinguish the vortex solution of the Abelian-Higgs model from the vortex solutions 
of our multi- Abelian-Higgs models, by using the word "normal" . 
In the Abelian-Higgs model, the Lagrangian density is 

c = -\f^f, v - |zvi 2 - \g 2 (- 2 - f 2 Y , (ci) 

where = d il A v — d u A^ is a field strength of the Abelian gauge field A^, a is a complex 
scalar field and / is its vacuum expectation value (a) = f. D^a is the covariant derivative 
of the scalar field 

Dpo = <V + igA^a, (C2) 

where g is the gauge coupling constant to the scalar field a. 

To obtain the classical solution in this theory, we impose the static, axially-symmetric 
ansatz: 

A = e v P(r) , (C3) 
a = p(r)e tnv , (C4) 

where the integer n is the winding number. We used the circular cylindrical coordinates r, 
ip, and z. 

We use the scale conversion x = gfr, P = gP and p = p/f. Therefore the energy density 



of per unit length of the z axis becomes 



£ = 2vr/ 2 / dx x 



POO 



i ^ , + ^-iV + ^ + ^V--^-i)-^"" 



, (05) 



2 \ x J \ x J x x 

where the prime (') denotes the derivative with respect to x. Asymptotic values are as 
follows: P(0) = 0, P(oo) = n, p(0) = and p(oo) = 1. We can write the following 
inequality for the energy 

roc 

£>2nnf 2 {p 2 )' dx = 2nnf. (C6) 
Jo 

This lower bound on the energy is the Bogomolnyi bound and it is saturated when p and P 
satisfy the following equations 

^ = -(p 2 -1), (C7) 

t = (08) 
p X 

These equations are the Bogomolnyi equations. 
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APPENDIX D: ACTION AND EQUATION OF MOTION WITH VORTEX 
ANSATZ 



In this Appendix, we show the details about the Bogomolnyi equations for the vortex 
configuration. We take the axially symmetric ansatz: 



Then we find 



a e = p e (r)e m ^, A v = P v (r) . 



V r a e = p' e e m ^ , V v a e = i{n e + (gP) t(e) - {gP) o{e) )p e e m ^ , 



(Dl) 



(D2) 



where the prime denotes ^, the derivative with respect to r, and (gP) v = g v Pv Thus the 
kinetic term of the scalar reads 



IT, |2 / /X2, ( n e + (gP)t(e) ~ {gP)o{e)f 2 
\Di°e\ = KPe) H 5 



while the Maxwell term becomes 



1 1(P') 2 

4 v ij ~ 2 r 2 



(D3) 



(D4) 



The total action can be rewritten as 



= 2tt f 
Jo 



drr 



+ \ E^-Ce)(^G 2 EW(^-CeO 



e,e'£E 



(D5) 



and this is no other than the energy density per unit length in the present static case. 
Varying this, we obtain the following equations of motion: 

— = ((gTGP) r"' )2 p« + £ - M , (D6, 



e,e'<=E 



7 e£E 



((E T GP) e -n e ) 2 T 

—o Pe\^ 



(D7) 



These second-order simultaneous equations can be reduced to the first-order Bogomolnyi 
equations: 

{E T GP) e -n e 



Pe = T 
P' 



Tj2(Pe-(e)(E T G) ev . 



(D8) 
(D9) 



e£E 
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APPENDIX E: ASYMPTOTIC PROFILE OF THE VORTEX 



We investigate the asymptotic behavior of the solution of (ID8|D9j) in this Appendix. To 



this purpose, first we introduce new variables p v (r) and R e (r): 

P v (r) =a v - p v (r) , p e = f e - R e (r) , (El) 
where the constant a v satisfies 

n e = (E T Ga) e . (E2) 

Next we prepare p-dimensional eigenvectors x^ (a = 1, • • • ,p — 1) for the (mass) 2 mass 
matrix for vector fields satisfying 

2(GEF)(GEF) T x {a) = (m (a) )V a) , for nonzero modes (E3) 

and g-dimensional eigenvectors for the (mass) 2 mass matrix for scalar fields satisfying 

2(GEF) T (GEF)X {a) = (m {a) ) 2 X {a) . (E4) 

Hereafter we restrict ourselves on the case with tree graphs treated in the text. Thus 
q = p — 1. The zero mode satisfies 

2(GEF)(GEF) T x w =0. (E5) 

The relations of two sets of eigenvectors are 

xV = ^(GEF) T xV, xM = ^jLGEFXV, (a ? 0) (E6) 

and we adopt the normalization convention: 

x (a)T x (a) = X (a)T x (a) = j (E?) 

Using the eigensystems, we can expand the variables by eigenvectors as 

Pv(r) = 5> (a) 4 a) , Re(r) = R(a)X e ] . (E8) 

(a) (a) 

Noticing R e (oo) = and p v (oo) = 0, the equations of motion ( 1D6|D7I) becomes at the 
asymptotic region, r — > oo, 

R (a)» + \ R {a)l _ ( m («))2 jR (a) = g ^ ^ 
24 



p(«)" _ IpW _ ( m ( a ))V a ) = , (E10) 
r 

and the Bogomolnyi equations ( 1D8|ID9I) become at the asymptotic region, r — > oo, 



1 ™, (a) 

Rio)' = («) , (Ell) 

— = -V2m {a) R {a) . (E12) 
r 

The solution of the above equations is 

R {a) = CK (m^r) , p (a) = V2C rK^m^r) . (E13) 

This result can be derived by using the following formulas for the modified Bessel function 
of the second type, such as K (z) and K\(z)\ 

K'^z) + -K' (z) - K (z) = , K'({z) + -K[(z) + K^z) = , (E14) 

{zK x (z))" - \zK x {z))' - {zK x {z)) = , (E15) 

z 

K' (z) = -K^z) , {zK x {z))' = -zK (z) , (E16) 

where the prime (') means the derivative with respect to z. 

More rough estimation can be done with the exponential function because 

K M « \hr e ~ z > for lar s e z ■ ( E17 ) 

V 2z 
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